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Abstract 

We consider Ising models in two and three dimensions, with short 
range ferromagnetic and long range, power-law decaying, antiferro¬ 
magnetic interactions. We let J be the ratio between the strength of 
the ferromagnetic to antiferromagnetic interactions. The competition 
between these two kinds of interactions induces the system to form 
domains of minus spins in a background of plus spins, or vice versa. If 
the decay exponent p of the long range interaction is larger than d-l-1, 
with d the space dimension, this happens for all values of J smaller 
than a critical value Jdp), beyond which the ground state is homoge¬ 
neous. In this paper, we give a characterization of the infinite volume 
ground states of the system, for p > 2d and J in a left neighborhood 
of Jc{p)- In particular, we prove that the quasi-one-dimensional states 
consisting of infinite stripes {d = 2) or slabs (d = 3), all of the same 
optimal width and orientation, and alternating magnetization, are in¬ 
finite volume ground states. Our proof is based on localization bounds 
combined with reflection positivity. 


1 Introduction and main results 

The problem of proving the emergence of periodic patterns in systems with 
competing interactions is ubiquitous in several areas of physics, biology and 
material science |2] , ranging from superconductor physics |T2] , micro-magne¬ 
tism ra, polymer suspensions [23], martensitic phase transitions [23], quan¬ 
tum Hall systems [H], to metal-oxide-semiconductor held-effect transistors 
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EH, nuclear matter [31] , and many others UU- In all these systems, experi¬ 
ments or simulations show evidence for the formation of remarkable patterns 
in suitable regions of the phase diagram. Examples are stripes, bubbles, zig¬ 
zag patterns, and columnar phases. The fundamental understanding of these 
phenomena is still in a primitive stage, mostly based on variational compu¬ 
tations and special assumptions on the structure of the low-energy states. 
There are just a few special cases where the periodicity of the ground state 
can be proved from first principles [SI 121 ESI ESI HO]. 

A particularly interesting and poorly understood phenomenon is that of 
periodic stripe formation [21 EHl ESI EE]- In a series of papers, this phe¬ 
nomenon was studied in Ising and related models with short range attractive 
and long range repulsive interactions. The method of block reflection posi¬ 
tivity led to rigorous proof the existence of periodic striped states both in one 
dimension [ISIEIIIE] and in certain two-dimensional models, including a toy 
model model for martensitic phase transitions [21] and a model of in-plane 
spins with discrete orientations and dipolar interactions [16]. However, the 
physically interesting case of out-of-plane spins with dipolar, or dipolar-like, 
interactions, which is of great importance for the physics of thin magnetic 
films, eluded any rigorous treatment so far. 

In a recent work 121 . we succeeded in computing the specific ground state 
energy of such a system, with power-law decaying repulsive interactions and 
decay exponent p > 2d in d = 2,3 space dimensions, asymptotically as the 
ferromagnetic transition line is approached. Our estimates allowed us to 
prove emergence of periodic stripe order in a suitable asymptotic sense, but 
they were not strong enough to fully control the ground state structure, or 
to prove breaking of rotational symmetry in the ground state. In this paper 
we extended the ideas of [ISIEQ] and prove that periodic striped states of 
optimal width are exact infinite volume ground states. Moreover we give 
a characterization of infinite volume ground states that are invariant under 
translations by one (for d = 2) or two (for d = 3) independent fixed lattice 
vectors. 

The setting is the following: consider Ising models defined by the formal 
Hamiltonian 


H{a) = -J - 1 ) 

(x,y> 


E 

{x,y} 


[a^ay 


1 ) 


( 1 . 1 ) 


where a G {±1}^ , d > 2. The hrst sum in (1.1) ranges over nearest neigh¬ 
bor pairs in while the second over pairs of distinct sites in For 
different values of the exponent p, this model is used to describe the effects 
of frustration induced in magnetic films by the presence of dipolar inter- 
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actions (p = 3) or in charged systems by the presence of an nnscreened 
Conlomb interaction [p = 1), as well as many other frnstrated systems 

[IlEllllElE10Elini[IIl[22lEZlEHlEniE2lE3lEaESlEHlllI]. 

In this paper, we choose the exponent p to satisfy the constraint p > 2d. 
As discnssed in a previons work [20], if J > Jc, with 


Jc-.= 


E 

yi>o, 


Vi 


( 1 . 2 ) 


then there are exactly two gronnd states, dx = +1 Vx G Z'^, and cXx = — 1 
Vx G Z'^. For J < Jc, the gronnd state is not nniform, and for J close to 
Jc it was conjectnred to be a periodic striped conhguration, i.e., a qnasi- 
one-dimensional periodic conhgnration of the form = sign(sin(7r(a:i + 

l/2)/h)), or translations or rotations thereof, for a snitable stripe width 
h G N. In this paper, we prove this conjectnre, and characterize the set of 
inhnite volnme gronnd states invariant nnder translations generated by d — 1 
lattice vectors. 

For simplicity, we restrict the discnssion to d = 2 and p > 4 from now on. 
Similar considerations are valid in d = 3 and p > 6 (or, in fact, for any d > 2 
with p > 2d). In Appendix]^ we explain how to adapt the proof to dimension 
three and higher. Let es(h) be the energy per site of compnted via ( |1.1[ ). 
We let h* = argmin;jgj^es(h), which is nniqnely dehned for almost alQ choices 
of J. We denote by qj = G {±1}^ , and we call it an optimal periodic 
striped configuration. Other Ah* — 1 optimal periodic striped conhgnrations 
are obtained from q* via translations and rotations. 

In order to state onr main resnlt, we also need to introdnce the following 
notions: the conhgnration s G {±1}^^ is called an infinite volume ground 
state if it is energetically stable against compactly snpported pertnrbations, 
that is, for any finite X C ifi, 

Hx{qxk) - Hx{sx\s) > 0, 

where Sx is the restriction of s to X, and 

Hx{qxk) = -J 5^(axCTy-l) + 

(x,y>: 

X,yex 

— J (f^x-Sy — 

xex, yex'^: xex, y&X<^ ' 

|x-y|=l 

^There are exceptional values of J for which es{h) has two minimizers, h* and h* + 1. 


Wqx G {±1} 


X 


(1.3) 


E 

{x,y}: 

x,yeX 


0+ E 


(axCTy - 1) 


|x-yP 


(o'x’Sy 1) 
lx - y\p 


. (1.4) 
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We shall say that two inhnite volume ground states are equivalent, if they 
only differ on a finite set. The equivalence class of a given inhnite volume 
ground state is called a sector. A sector is trivial, if it contains only one 
element. In terms of these notions, onr main result can be snmmarized as 
follows. 

Theorem 1. There exists e > 0 such that, if — s < J < Jc, then the 
optimal periodic striped configurations are infinite volume ground states, and 
their sectors are trivial. 

This result is a corollary of a quantitative lower bonnd on the energy 
of spin conhgurations, which will be formnlated in Theorem below, after 
having introduced a few more dehnitions. Our quantitative bonnds also allow 
us to characterize the inhnite volnme gronnd states that are invariant nnder 
translation by a vector n G 

Theorem 2. Under the same conditions as Theorem^ any infinite volume 
ground state s that is invariant under translation by a vector n = (ni,n 2 ) G 
is characterized by the following property: there exists an “interface” of 
finite width, of the form Xfc^^fc 2 (n) = {x G : x ■ nj_ G [ki,k 2 \}, where 
n_L = {—n 2 ,ni) and ki < k 2 are two integers, such that s coincides with 
two of the optimal striped configurations on the two infinite components of 

(Xfci,fc2(n))''. 

Let us now introdnce a few more dehnitions, which are reqnired for the 
formulation of onr qnantitative lower bonnd on the energy of a generic spin 
conhguration. 

1.1 On good and bad 

1.1.1 Contours and corners 

Given q_ G {±1}^^, we let A = {x G Z^ : (a)x = —1}, and r(A) be its 
boundary, i.e., the nnion of bonds of the dnal lattice (Z^)* separating a point 
X G A from a point y G A'^. At every vertex of r(A) fl (Z^)*, there can be 
either 2 or 4 sides meeting. In the case of 4 sides, we deform the polygon 
slightly by “chopping oh” the edge from the squares containing a — spin; see 
Fig.[T} 

After the chopping, r(A) splits into disconnected polygons ri,r 2 ,..., 
which are called contours. The points where two orthogonal portions of 
a contonr meet are called corners. The sites of the dnal lattice where a 
non-trivial chopping operation took place correspond to two corners. We 
denote by Nc{Ti) the number of corners of Tj, and Nc{A) = AfdTi). Note 
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Figure 1: In the case that 4 sides of the closed polygon r(A) meet at a 
vertex v, we slightly deform r(A) so that the two squares containing a — 
spin become disconnected from the vertex itself. Case (a) represents the 
situation where the minus spins are located at NE and SW of v, before and 
after the “chopping”. Case (b) represents the situation where the minus spins 
are located at NW and SE of v, before and after the “chopping”. 


that, if a is a compactly supported perturbation of a*, that is, if the set 
{x G : (a)x 7^ (a*)x} is finite, then Nc{A) < +00. 

1.1.2 Tiles 

Given an integer i, we pave with tiles of side £ and, given a tile T, we 
denote by Ft (A) the union of the bonds in F(A) that separate a point in AflT 
from a point in A'^. Note that this convention assigns uniquely every bond in 
F(A) to one of the tiles. The connected components of F 7 ’(A) are contained in 
the contours Fi,..., and are denoted by F^ j, i = 1,..., r^. Given a maximal 
straight portion of F^,*, we assign to each of its two endpoints a “number 
of corners”, which can be either 1/2 or 0, depending on whether or not the 
given endpoint coincides with one of the corners in Uj>iFj. This assignment 
induces a notion of “number of corners in the tile T”, to be denoted by Uc(T), 
which is the sum of the number of corners of all the endpoints of the straight 
portions of F^^j, with z = 1,..., tt- Note that nc{T) can be either integer or 
half-integer, and J2t'^c{T) = Nc{A). 

1.1.3 Bad tiles and good regions 

We now identify the tiles Ti such that either nc{Ti) > 0, or they contain a 
square of side £/5 completely contained in A or in A'^, to be called hole. We 
call these tiles bad^ and we let Afs be their number. For later convenience, 
we also let be the number of bad tiles containing a hole. The connected 
components of the complement of are denoted by G*, z = 1,..., A/g, 

and are called the good regions. By construction, any of these connected 

^We shall choose £ large compared to the optimal stripe width h*, which explains why 
we expect a hole to be energetically unfavorable, hence bad. 
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components contains portions of contours in r(A) that are all straight with 
the same orientation, and have no corners. We denote by Fcj the union of 
contours in r(A) contained in Gj. If the elements of Fc; are all vertical, 
we will say that Gi is vertically striped, and horizontally striped otherwise. 
Consider a good region Gi that is vertically (resp. horizontally) striped. We 
say that i? is a “rectangular portion of stripe” in Gi, if i? is a rectangle 
completely contained in A or in A'^, with its two vertical (resp. horizontal) 
boundaries both belonging to Fg.. We also define the distance between its 
two vertical (resp. horizontal) boundaries to be the width of R. Finally, we let 
Ah{Gi) denote the area of the union of all the rectangular portions of stripes 
of width h in Gi. Note that, if u is a compactly supported perturbation of 
q*, then ^h{Gi) is hnite, \/h ^ h*. 

We are now in the position of stating our quantitative lower bound on 
the energy Hx{sx\q^) of a spin configuration s with a* boundary conditions. 


Theorem 3. There exist positive constants Go, Gi, e such that, if J^ — e < 
J < Jc and Goh* < £ < {Go{Jc — then for every s G {±1}^^ and every 

finite set X (ZlA, 

Hx{sx\cL*) > Hx{q*x\q*) + Gr{N, + {R-jYt£^H^Ml°^^) 

Mg 

+ ~ ^s{h*))Ah{Gi), (1.5) 

h^h* i=l 

where Nc, and Gi are, respectively, the number of corners, the number 

of bad tiles containing a hole, and the good regions, associated with the infinite 
spin configuration a = {sxyT*x‘=) coinciding with Sx on X and with a* on 
X^, defined via tiling with sguares of side length £ as described above. 


Remark. Since h* ~ {Jc ~ J) ^ {Jc — J) ^ for small Jc — J 


(compare with Eq. (2.30) below), the condition on £ can be fulfilled for J 
close to Jc. Note also that es{h*) ~ {Jc — J)(p“2)/(p- 3) small Jc — J, which 
agrees with the factor multiplying in the second term on the right side 


of (1.5) 


Remark. The prefactor 1/2 in the second line of (1.5) can be replaced 
by any number less than 1, at the expense of modifying the constants Go and 

e. 


Theorem implies, in particular, that of is an infinite volume ground 
state, and that every state a that is a compactly supported perturbation 
of it, is not a ground state, simply because any such state necessarily has 
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corners and, therefore, by ( |1.5 ), it has strictly larger energy than a*. This 
immediately implies Theorem 

In order to see that also Theorem is a conseqnence of Theorem 
consider an infinite volnme gronnd state s that is invariant nnder translations 
by an integer vector n. Let A C be a sqnare box of side L, and note that 
the energy price for changing the bonndary conditions from s to a* scales 
like the bonndary, that is 


\HA{a^\s)-H^{a^\a*)\<2{J + J,)\dA\ ( 1 . 6 ) 

for any Using this ineqnality and the very definition of infinite volnme 
gronnd state, we have 


HA{al\a*) + 2{J+J,)\dA\>HA{al\s) > (1.7) 

> HaUj^Is) > HAisj,\a*) - 2{J + Je)|M|. 

Now we apply Theorem]^ thus obtaining 

C'i(iV,+(J,-J)l^£W^°''=)+- 5^(es(h)-e,(h*))A,(G,) <4(J+Jc)|aA|, 

h^h* i=l 

(1.8) 

where and Gi refer to the configuration (s^,a)^c). In particular, 

Nc is bounded by (const.)L. Since every corner not at the boundary of A is 
repeated with period n, there can be at most a finite number of them (modulo 
translations by n) independently of L. This means that these corners are all 
contained in a finite strip Xfc^fc 2 (n), as claimed in Theorem]^ with fci,/c 2 
independent of L. Similarly, we can argue that the holes and the stripes 
of width different from h* are all contained in a hnite strip Xki^k 2 {'<^)- This 
concludes the proof of Theorem]^ in the case that h* is unique. As observed 
above, there are exceptional values of J for which es(h) has two minimizers, 
h* and h* + 1. In these cases, the discussion above leaves open the possibility 
that on one of the connected components of the stripes are not 

all of the same width. However, this cannot be the case: by proceeding as 
in m Section III.D], one can prove that each pair of neighboring stripes of 
widths h*,h* + 1 gives an extra positive contribution to the energy per unit 
stripe length. Therefore, pairs of stripes of different widths are all contained 
in a finite strip X^^ Theorem follows. 

The rest of the paper is devoted to the proof of Theorem 
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2 Proof of Theorem |3| 


The proof of Theoremis divided into several steps, and uses many notations 
and ideas introduced in I2D1, which will be recalled here. As a preliminary 
step, we reduce to plus boundary conditions, which allows us to use the 
droplet formulation for the energy as in [20], see (2.3) below. Once the 


energy is expressed in terms of droplets, we can localize the energy in the 
bad tiles and good regions, by proceeding in a way analogous to [201. The 
key technical novelty of this paper, as compared to 121 . is an efficient way 
of estimating the energy in the good regions, which may have a complicated 
geometrical shape. The crucial estimate is summarized in Lemma below. 


whose proof is given in Section 2.4 


Before entering the proof, we recall some notation: we let r = 2(J — Jc), 
which is assumed to be negative and small, in absolute value. We recall that 
h* is of the order and the specific energy of an optimal striped 

configuration, es{h*), is negative and of the order (compare with 

( 2.30| ) below). 


2.1 Reduction to plus boundary conditions 

The a* boundary conditions, while natural in the perspective of proving 
uniqueness of the ground state, are not particularly convenient for using the 
droplet representation of 121 . However, a few simple algebraic manipula¬ 
tions allow us to reduce to the same boundary conditions of [20] (that is, 
plus boundary conditions) in a suitable enlarged box. To see this, rewrite 
Hx(^x\S‘) as 

Hx('Lx\'L') = Hx(s:x\'f) + ( 2 - 1 ) 

where A is a square box of side L, which we choose to be divisible by 2h*, 
configuration on A whose restriction to X C A (resp. A \ 
X) coincides with ax (resp. and iL^®’^(a^) = H/^{aj^\a^^) is the 

Hamiltonian with periodic, rather than a*, boundary conditions (here 
is the periodic extension of aj^ over Z^). 

Now, = es(h*)|A| and Hj^^\ax,Qi*j^\x) can be further rewritten 

in terms of a Hamiltonian with plus boundary conditions: 

Hl'^'iax, q:X\x) = Jp^Am Aa\x)^ ) (2-2) 

where Am is a square box of side LM, to be thought of as the union of 
copies of A, and ^ symbol for the configuration on Am 







obtained by juxtaposing copies of one in each of the copies 

of A. Moreover, H^{ax) = Hx{oix\QL^) indicates the Hamiltonian with 

plus boundary conditions (here a"*" is the uniform inhnite spin conhguration 

consisting of plus spins everywhere). Finally, for later reference, we introduce 

the shorthand for the spin conhguration {o^x, and u for 

the inhnite one coinciding with on Am and with on the complement. 

From now on we shall consider Ma and u hxed once and for all. 

—— 


2.2 Localization 


We now re-express in terms of the droplets representation intro¬ 

duced in 


Using the notation introduced in Section 1.1.1 


we let A be the 

region of minus spins associated with u, F(A) its boundary, and Fi,..., F^ 
the corresponding contours. We also denote by 0(A) the collection of con¬ 
tours, 0(A) = {Fi,...,Fr}. Note that A, F(A), 0(A) and r are hnite, 
because of the plus boundary conditions. As in PDl, we denote by Si the 
maximal connected components of A, and by T>{A) their collection. Given 
S G T>{A), we also let F(5) be the boundary of S, and Nc{S) its number of 
corners. In terms of these notations, we can re-express the energy of as 


= H iri 


re©(A) 






(2.3) 


S&V{A) 


S,5'e:V(A) 
Sj^S' 


where 

U(S):=-2J2 E 


xe<5 y£l,2\S 


|x-y| 


W(i,S') :=4EE 


xGS yGS' 


x-y| 


(2.4) 


Let us consider the partition P of Am de hned by the bad tiles and good 
regions of Am, in the sense of Section [1.1.3[ V = U We now 


localize the energy in the elements of V, by proceeding as in [2D1 Section 
3]. More precisely, we derive a lower bound on the energy in 

the form of a sum of local energies Eq{Bq), each depending only on the 
“bubble conhguration” Bq within the region Q q V. The notion of bubble 
conhguration was introduced in [201 Section 3] and is recalled here: given 
Q E V and S E V, we denote by Fq( 5) the portion of F(5) belonging to 
Q. Moreover, if (5 q = 6 f] Q, we dehne Sq\ ..., to be the maximal 

connected components of Sq, and ..., portions of Fq((5) 


belonging to the boundary of Sq \ ..., Sq"^''‘'’\ respectively. We shall refer to 


Q 


the pair (5 q\Fq^) as to a bubble in Q originating from S. We shall indicate 
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Figure 2: On the left: a square region Q (in light grey, with dotted boundary) 
and the droplets having non-zero intersection with it. On the right: the 
corresponding bubble conhguration Bq after localization in Q. Note that a 
single droplet can give rise after localization to more than one bubble: e.g., [82 
and ( 8 ^ originate from the same droplet, and similarly for ( 8 ^ and ( 84 ^. Note also 
that the contour of a bubble does not necessarily coincide with the boundary 
of its droplet: in general, it is contained in it, and may even be disconnected 
(as in the case of ( 8 ^, whose contour consists of two disconnected portions). 


by Bq{6 ) the set of bubbles in Q originating from 6 , and by Bq = Us£vBq{6 ) 
the total set of bubbles in Q. (See Fig. i) 

Given a bubble [8 = ( 5 / 3 , F/ 3 ) G Bq, we also dehne its localized self-energy 
as 

“«(//) = - E E (2-5) 

where is the distance between b and the bond 6 ' G F /3 facing it^in (8, 

^ The notion of “bond facing b in /3” is defined as follows. Let us suppose for definiteness 
that 5 G T/j is vertical and that it separates a point xt G on its immediate right from 
a point = X(, — (1,0) ^ on its immediate left. We say that h faces b' G T^, and vice 
versa, if: (i) 6' is vertical; (ii) b' separates a point Xf,/ G Sp on its immediate left from a 
point yb' = Xf,' -|- (1,0) ^ Sp on its immediate right; (hi) the points Xf, and Xf,' are at the 
same height, i.e., [xf ,]2 = [xt,/] 2 , and all the points on the same row between them belong 
to 6 p: in other words, xt + (j, 0) G Sp, for all j = 0,..., [xf,']i — [xf,]i. An analogous 
definition is valid for horizontal bonds. If b does not face any b' GTp, we say that b faces 
the boundary of Q. 


10 


















































































if there is such a bond, and is inhnite, if b faces the boundary of Q. In terms 
of these notions, our localization bound takes the following form: 

A/fl A/g 

> E + x; EcJBc.), (2.6) 


where 

et.(Bt.) = 5 ^ i 2 j|ni+ UT,m + ^ E v)+ 2 ‘-'’/V(t.), (2.7) 

/36Bt, /3./3'eBT,, 

Eg,(Bg,) = Y1 + “a(/^)] + § E V). (2.8) 

/366Gi P,fi'eBG,, 

/3^/3' 


and the * on the last sum indicates the following constraint: if the bubbles in 
Gi are all vertical (resp. horizontal) we only sum over pairs of bubbles that 
do not overlap after arbitrary translations in the vertical (resp. horizontal) 
direction. 


Proof of (2.6). We start from (2.3). The goal is to bound it from below 
by a sum of terms, each of which is localized in an element Q of V. The 


hrst term on the right side of (2.3) is already local, i.e., it can be rewritten 
exactly as ’''^^ich leads to the corresponding terms in 

(2.7) and (2.8). The interaction W{S,6') on the right side of (2.3) can be 


rewritten as 


»'(^A')= E E 'r(i^,V), 


Q,Q'&P P&Bq{&) 
/3'6Bq/(5') 


which is bounded from below by dropping the terms with Q ^ Q' (recall that 
the interaction is positive), so that 


1 E E 'r(ij.V). 


(2.9) 


d,5'&V{A) 

S^S' 


Q&P p,p'eBQ 


and the ** on the sum indicates the constraint that the two droplets 5, 5' in 
P(A), which 5^, 5pi belong to, are different from each other, 5 7 ^ 5'. 

Regarding the second term on the right side of (2.3), we bound it from 
below by using [201 Eq.(2.9)], that is 

u(S) > - E E + 2 ‘-'iv,('i)+4 E ^ 


ber((5) n^o 


{x,y}eP((5) 


x-y|P 


( 2 . 10 ) 
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where db{ 6 ) is the distance between b and the bond b' facing it in 6 , and V{ 6 ) 
is the set of nnordered pairs of distinct sites in 6 snch that both and 

cross at least two bonds of r((5). Here is the path on the lattice 
that goes from x to y consisting of two segments, the first horizontal and the 
second vertical. Similarly, is the path on the lattice that goes from x 

to y consisting of two segments, the first vertical and the second horizontal 
(note that the two paths can coincide, in the case that x* = i/i for some 
i€{l,2}). 


The first term on the right side of (2.10) can be bounded from below as 


E E 

ber((5) n^o 


min{|ni|,db(5)} 


n P 


> E E “o(/3). (2.11) 

Q&ViA) l3eBQ{S) 


which, after summation over 5, leads to the terms in (2.7) 


and (2.8). The second term on the right side of (2.10) is local and, after 


summation over 5, can be rewritten as 2^ 2 ndTi) (recall that the good 
regions have no corners), which leads to the last term on the right side of 


(2.7). Moreover, the sum over 5 of the last term on the right side of (2.10) 


can be bounded from below by an expression similar to the right side of (2.9), 
namely 


E^ E 

S&V(A) {x,y}e(P(<5) 


|x-y| 


1 

> - 
- 2 


t 


5 ; 5 ; 


( 2 . 12 ) 


Q&V 13,0'&Bq 


where the f on the sum indicates the constraint that 60 , 60 ' belong to the 
same droplet 6 G T>{A), and all the pairs of points (x, y) in 60 x 60 ' are such 
that {x, y} G V{ 6 ). Combining the right sides of (2.9|) and (2.12) we obtain 


tt 


5 ; ^ W(Sf,Sr,), 


Q&V 0,0'&Bq 

where the ff on the sum indicates the constraint that: either 60 , 60 ' belong 
to different droplets in V{A), or, if they belong to the same droplet in 6 , 
they are such that all pairs of points (x, y) in 60 x 60 / are in V{ 6 ). Finally, 
note that: if Q is a bad tile (which is a convex region), then the constraint ff 
is automatically realized (i.e., it can be dropped), which leads to the second 


term on the right side of (2.7); if Q is a good region, then the constraint ff 


is easily seen to be weaker than the one indicated by * in (2.8), which leads 
to the last term on the right side of (2.8). □ 
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2.3 Lower bounds on the localized energies 


In this subsection we state the key lower bounds on the localized energies in 
the bad and good regions, and prove that they imply Theorem Recall the 
definitions of Mb-, and Ah{Gi) given in Section 1.1.3[ Recall also that ^ 

is the side length of the tiles, which enters the definition of the partition V. 


Lemma 1. There exist positive constants Cq, Ci and e such that, if —e < r < 
0 and E > c^h*, then the energy Eg of any good region G eV satisfies 

Eg{Bg) > e,{h*)\G\ - c,\T\\dG\ + ^ (e,(h) - es(h*))R,(G'). (2.13) 

^ h^h* 

Lemma 2. There exist positive constants cq, C 2 and e such that, if —e < 
r < 0 and Coh* < i < (co|r|)“^, then the energy Ex of any bad tile T E V 
satisfies 

Et{Bt) > f e,(h*) + C 2 [n,{T) + , (2.14) 

where Xhoie(7") is egual to 1 if T contains a hole, and 0 otherwise. 


Lemma [T] is one of the main technical novelties of this paper, and its proof 
is described in detail in the next section. The proof of Lemma is simpler, it 
is an extension of the bounds worked out in [20], and its proof is postponed 
to Section [231 

Combining the two lemmas, we can easily derive Theorem]^ In fact, note 
that every portion of the boundary dGi of a good region Gi eV is adjacent 
to a bad tile, so that Ym=i — 4£A/b. Using this observation, together 
with the fact that every bad tile either has a positive number of corners (at 
least 1/2) or a hole, and plugging (2.13) and (2.14) into (2.6), we obtain 

> es(h*)|AM| + I (iV,(A) + 

e,{h*))A^,{Gi) 

— 4ci|r|£j. 


Ma 

h^h* i=l 


(2.15) 


+Mb[-^ 


C2 


mm < - 


Now, if 


8ci 

C2 


|-l/(p-3) 




C2 

16ci 


T 


1-1 


then the expression in parentheses that multiplies Mb is positive, and we can 
drop it for the purpose of a lower bound. Via the use of (2.1)-(2.2), we thus 
arrive at Eq. (1.5). 
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Figure 3: A good region with its bubble configuration. The regions in dark 
grey are droplets, while those in light grey are their complement. Note that 
not all of the droplets are rectangular, due to boundary effects. E.g., the 
droplets 5i and 62 are not rectangular. The arrows indicate the direction in 
which the vertical boundary segments move under the deformation described 
in the text. 


2.4 Proof of Lemma [U 

2.4.1 Deforming the good regions 


We now discuss the proof of Lemma Let us consider a good region G &V, 
and let us assume without loss of generality that it only contains vertical 
stripes. We first want to slightly deform the domain G and correspondingly 
change the bubble configuration within, in order to make all the bubbles 
rectangular: here we call “rectangular” a bubble {3 = ( 5 / 3 , T^j) such that 5^ 
is a rectangle and r /3 is the union of its two vertical sides. Note that, in 
general, not all the bubbles in G are rectangular, due to boundary effects (in 
fact, the boundary can partially “cut” a portion of the rectangle, without 
disconnecting it, see e.g. the droplets 5i and 82 in Fig. [^. 

In order to describe the deformation of G, think of the vertical boundary 
of G as a union of segments Si of length £, induced by the tiling described in 
Sections |1.1.2| and |1.1.3[ By construction, every boundary segment Si faces 


a portion Si of length ^ of the boundary of a rectangular bubble in Bq that is 
closer than 2i/5 to Si itself. We now deform the boundary dG continuously. 
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Figure 4: The deformed good region corresponding to the good region of 
Fig-i after the moves of the vertical boundary segments indicated by the 
arrows. Note that after the deformation all the droplets are rectangular, 
and all the connected components of the complement of the droplets (the 
connected light grey regions) are rectangular as well. 

by moving the segments Si towards the interior of G, in such a way that they 
coincide with Sf, see Fig. 

In this way we increase the boundary of G by at most \dG\ itself, and we 
increase the energy by at most |r| ■ |9G|, which is acceptable for our purposes. 
We shall denote by G' the new region obtained from G by the deformations 
we just described, and the new bubble conhguration Be- We have 

Eg{13g) > Eg\Bg') - |r| ■ I^GI, (2.16) 


as just argued. 

2.4.2 Slicing and bounding the energy of the good region 

The deformed region G' obtained in the previous step is a union of connected 
horizontal slices gj of height £, with j = 1,..., iV^, as shown in Fig. To 
each slice we associate a sequence of integers (hi, tci,..., tCn-i, h^), where n 
is the number of bubbles in Bg' intersecting the slice, hi, ^ 2 ,..., h„ are their 
widths, ordered from left to right, and Wi,W 2 ,... are the spacings between 
the hrst and second bubbles, second and third, etc. 
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Figure 5: The six slices corresponding to the deformed good region of Fig. 
For ease of visualization, the slices are drawn slightly detached from each 
other. As an example, in the hrst slice we attached the labels hi and Wi 
indicating the widths of the stripes (dark grey regions) and of the rectangular 
regions separating two subsequent stripes from each other (light grey regions). 
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We also denote by Si, with i = 1,..., W, the maximal connected seg¬ 
ments in the intersection between the horizontal boundary of G' and the 
boundaries of the droplets in Bq' (note that the boundary of a rectangular 
bubble does not coincide with its contour: rather, it consists of four seg¬ 
ments, two horizontal and two vertical). Note that segments come in pairs: 
one can say that two segments form a pair if they belong to the boundary of 
the same droplet. Moreover, to each segment Sj we associate its length h{sj) 
and two spacings Wi{sj), W 2 {sj), which are the horizontal distances from the 
next droplets (“next” by following the boundary of G') to the left and to the 
right of Sj. Note that if Sj touches a corner of dG', say on its right side, then 
there may not be any droplet to its right: by following the boundary one 
may hnd that the next segment to the right could actually have the same 
horizontal coordinates, in which case we will assign the value -foo to W 2 {sj) 
(and similarly for Wi{sj) in the case of the next segment on the left); see 
Fig. 

The key ingredient for the proof of Lemma is the following. 

Lemma 3. Given G' and Be as above, we have: 


Na 


Eg'{Bg') > ^ eoo(h{ 


U) 


w; 


i=i 


N, 

’ h{Sj),W2{Sj)), (2.17) 

i=i 


where h^i\w[^\ ... ,hn} is the sequence of widths and spacings associated to 
the slice Qj, and the functions and f are defined as follows: 


^ 00 (^ 1 ? ^ 1 ? • • • 7 hjfj 


*=i xez2\{o} 

+1 E W'ft.’Cft)). 


min{|xi|, hi} 


xF 


+ 




(2.18) 


where U = {(a;,0) E T? : 0 < a; — X]I=i(^fe + '^k) < K} and C{li) is the 
smallest infinite vertical strip containing U, that is C{li) = {(x, y) Elfi 0 < 
^ ~ + '^k) < hi, ye Z}; moreover, 

f{wi,h,W2) = ( 2 . 19 ) 

where Cf = {{x,y) e : —h < x < 0,?/ > 0} and Qiui,h,w 2 = ^ • 

X < —wi — h OT X > W 2 , 2/ < 0}. 
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Figure 6: The red segments Sj (color online) for the bubble configuration 
of Fig. 1^ in the deformed good region G'. The segments can be naturally 
grouped in pairs: two segments form a pair if they belong to the boundary 
of the same droplet. E.g., si and sn are paired, sn and S 20 are paired, 
etc. Every segment Sj comes with its length h{sj), and with two spacings, 
the left spacing Wi{sj) and the right spacing W 2 {sj), corresponding to the 
distances to the closest droplets to its left and to its right, by moving along 
the boundary. If one of the endpoints of sj is a corner of the boundary of 
G', then Sj may not have any droplet to its left or right, in which case we 
let the corresponding spacing to be inhnite. In the example in the hgure, 
rai(si) = tni(si 7 ) = W 2 {ss) = W 2 {sq) = + 00 , and all the other spacings are 
hnite. 
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Figure 7: A pictorial representation of the droplets involved in the definitions 
of f{wi,h,W 2 ), which represents the interaction energy of the infinite half 
strip of width h, with the two quarter spaces to its lower left and lower 
right. The spacings wi and W 2 represent the horizontal distances between the 
inhnite half strip and the quarter space to its left and to its right, respectively. 


Remark, e^o is the energy per unit vertical length of an infinite vertically 
striped conhguration, and f{wi, h,W 2 ) is the interaction energy between the 
droplets in Fig. 

Proof of Lemma First of all, note that the contribution 

5; i2j|ni+«G'(/3)] 


to the energy Eg'{Bg') (see (2.8)) is identical to the corresponding contribu¬ 
tion in \ ■ ■ ■ ^hn])-, i.e., to the one arising from the hrst 

of (2.18). Therefore, all we have to prove 
interaction terms in Eg'{Bg') and in 


two terms on the right side 
is that the difference between the 


^ ^oo{hi\wi \ ■ ■ ■, hn]) is bounded from below by the last sum in (2.17). 

Let us focus on a given slice gj and on the intersection of a given bubble 
(3 with this slice. The interaction energy of this portion of bubble with all 
the other bubbles in G', as it appears in Eg'{Bg'), is 


I E ms^ngj.s^,), (2.20) 

/3'eBG, 


where we recall that the * on the sum indicates the constraint that the 
bubbles (3' intersecting f3 after vertical translations should not be included in 
the sum. Eq. (2.20) can be bounded from below by summing only over the 
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bubbles /3' having non zero intersection with gj: 


J2.20| > - 


W{5png,,6p,). 


The term eoo{hi’ ,Wi\ ..., hn.) on the right side of (2.17) contains 

a term of that form, with the difference that <5/?/ is replaced by the inhnite 
vertical strip of the same width containing it. In fact, the interaction term 
in (2.18) satisfies 






N„ 


N„ 


'5^ = E w(Sf,ngi,c[h)). 


i=l 


i—1 P&Bqi P'&Bq! 


( 2 . 21 ) 

where is the analogue of the set U defined after (2.18), corresponding to the 
widths and spacings \ ..., hnj associated to the slice gj. Therefore, 

what remains to be proved is that 


Na 


Ns 


EE E W{ 6 pngj,C{S 0 ,)\ 6 p^) 

i=i 


j —1 PgBq/ P'^Bq/ 

P'ng.j^m, p'^p 


Let us rewrite the left side of (2.22) as 


E E E 


( 2 . 22 ) 

(2.23) 


P'&Bg> j=i,-,Ng p&Bg! 
5,-n/3V0 P^P' 


Note that the horizontal boundary of E(5/3') \ 5/3' consists of two segments 
Sj,Sfc with h{sj) = h{sk) and Sj above Sk (these are the pairs of segments 
mentioned in Section 2.4.2 see Fig. |^, and there is a one-to-one correspon¬ 
dence between summing over (3' and summing over these pairs of segments 
Sj , Sk (which is of course the same as summing over all segments). More¬ 
over, the set E(5/3') \ dpi is equal to the union of two sets, each of which is a 
translation, and one also a reflection, of that is 


E(5/3') \ 5/3' — U (2.24) 

for translations rj, and r|/, and r denoting reflection about the a:i-axis. 
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To conclude, it is enough to note that 


LJ LJ ^ 9i — '^p'Qw-i{sj),h(sj),W2isj) ^ 'l'i3'^Qwi{sk),h{^k),W2(sk)y ( 2 . 25 ) 

i=l,...,Ng I3 &Bqi 

9in/9V0 


which implies the desired inequality (2.22), and thus completes the proof of 
Lemma HI □ 


2.4.3 Reflection positivity 

We now show how to bound from below the right side of ( |2.17[ ), and how to 
use the resulting estimate to conclude the proof of Lemma [Tj The key step 
is to use reflection positivity to obtain a lower bound on the Coo term. This 
is an application of the block reflection positivity for one-dimensional spin 
systems worked out in [15]. The result is the following: 

n n—1 

eoo{hi,wi,...,hn) > r + ^hies(hi) ^'^Wie^{wi). (2.26) 

i=l i=\ 

To see this, note that eoo(hi, wi,..., hn) = lim^-^oo where: 

1. given a spin conhguration on = [1, L] fl Z, 

L 

^a“(^aJ = -J “ 1) + - j) , (2.27) 

i=l l<i<j<L 

with 

vl{x)= ^ {{x + nLf + , (2.28) 

and ai+i = crp 

2. the spins in the conhguration are equal to —1 on the intervals 

{x G Z : 0 < X — + Wk) < hi}, and -|-1 otherwise. 

Now, HYi^Lhi hn) ^ one-dimensional spin Hamiltonian with a reflection 
positive long-range interaction and periodic boundary conditions, of the class 
considered in m IT6] . Therefore, we can apply the chessboard estimate 
proved e.g. in the Appendix of [16]. As a result, using [m Eqs. (A4)- 
(A5)] and recalling the fact that the spin configuration consists of 

blocks of alternating sign, of size hi,wi,..., hn,Wn, with Wn = Wn{L) = 
L - {hi + wi -\ --h hn): we get 

n 

- '^{hies{hi) + WiCsiwi)) , (2.29) 

i=l 
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where es(h) is the energy per site (as computed from in the limit 

L —)■ cx)) of the inhnite periodic conhguration consisting of blocks all of the 
same size h, and of alternating sign. Note that es(h) is the same as the one 
defined in Section[^for the two-dimensional model. Finally, to go from (2.29) 
to ( 2.26D , observe that \m\L^ooWn{L)es{wn{L)) = r. This follows, e.g., from 
the explicit expression of es(h), derived in (JOl Appendix A]: 


“ ft + fc 


da 


(1 


o-a\2 


(1 _ tanh ^) = I + ^ + 0 (h 

(2.30) 

for large h, where /n^(a) is the inverse Laplace transform of the function Voo{,x) 
in (2.28), i.e., the function such that Voo{x) = da Wx > 0 , and 


Ap is a suitable constant. 


Remark. From Eq. (2.30) it follows straight away that the optimal 


stripe width is h* = {{p — 2 )Ap|r| _|_ o(1 )) as r —)■ 0 , and also that 

p-3 
p—2 h* 




2.4.4 Putting things together 


Plugging (2.26) into (2.17) gives 


EciBG') > \G'\e,{h*)+iTNg 

Ng rij Uj-l 

- e^{h*)) + - es{h*)) 

j=l k=l k=l 

Ns 

- Y '«^ 2 (sj)) . (2.31) 

1=1 


Now, £Ng < ^\dG'\ and it remains to show that the sum of the last two lines 
can be bounded from below by (const.)r|5G'|-|-^ —es{h*))Ah{G). 

From the definition of f{wi,h,W 2 ) it easily follows that it can be bounded 
independently of h, as 


f{Wi,h,W2) < Y 


G2 

0—4 
W- 

i=l,2 


(2.32) 


for a suitable constant 6 * 2 . Moreover, from (2.30) it follows that 


es(w) -es{h*) > + - 

VJP ^ w 


(2.33) 
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for all ty > 1 and a suitable > 0. Note that the left side is non-negative, 
while the right side may be negative. A simple consequence of (2.33) is that 

^ < \t\w + (C'3|r|"^)^w(es(w) - es{h*)). 


vjP- 


(2.34) 


Using this bound in (2.32) implies that the last line of (2.31) can be bounded 
as 

N, 

< (2.35) 


i=i 

C'a 

< — 
- 


Ns ^ 

+ iC 3 \T\~^)^Wi{sj){es{wi{sj)) - e^{h*)) 


j=l i=l,2 


where the prime on the sum indicates the constraint that Wi{sj) < oo. Now, 
every spacing Wi{sj) appears twice in the sum above (because every spac¬ 
ing is to the left or to the right of two different segments Sj , Sj /), hence 
<2|9G'|. Similarly, 


Ns 




Wi{sj){es{wi{sj)) - e^{h*)) < 2^ ^ - e^{h*)). 

j=l i=l,2 j=l k=l 

(2.36) 

Therefore, if £ > 


Eg'{Bg') > |G'|e,(h*) + T\dG '\(i + 2^) + 


(2.37) 


n-i — l 




Aj) 


(i)/ 




j=l k=l 


k=l 


To complete the proof of Lemma note that |G'| < |G|, and \dG'\ < 2|c?G|, 
as already argued above. □ 


Remark. The proof of (2.37) is valid for bubble conhgurations a bit 


more general than those considered here: in fact, we never used the fact that 


Bg' has no holes, in the sense explained in Section 1.1.3 The only property 


we really used is that the bubbles in Bg' are all rectangular with the same 
orientation. 


2.5 Proof of Lemma 

We proceed similarly to [201 Section 3]. The first step is to estimate the cost 
of erasing the bubbles with corners. Write ndT) = with z/c(/9) 
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the number of corners associated with the bubble /3, which may be an integer 
or a half-integer. Consider a bubble with UciP) > 0. Dropping the positive 
interaction of this bubble with the others, its contribution to the energy is 
bounded from below as 

2J|r^| + UTi^) + > r|r^| + 2^-"^u,{/3) . (2.38) 


Note that, in order for to be very long, the number of corners must be 
sufficiently large: in formulae (see [201 Eq.(3.10)] and following lines). 


\T0\<2i + 2iu,{P) . 


(2.39) 


If, as we are assuming, > 1/2, then -|- 1 < Sudd), so that I'dd) — 
|r^|/(6f'). Inserting this back into (2.38) gives 

2J\r0\+UT{d)+2^-^/^udd) > 2-P/^^^(l-6-2P/^\T\i^ +2-2udd) • (2.40) 


The first term on the right side is positive, and, therefore, can be dropped 
for a lower bound, if £ < (6 ■ 2^/‘^\t\)~^. Therefore, denoting by St the subset 
of Bt consisting of all the bubbles without corners. 


Et{Bt) ^ Et{St) + 2 "^ndT). 


(2.41) 


In order to estimate the energy of the corner-less conhguration St we 
proceed exactly as in the proof Lemma [1} Assume that the contours in St 
are vertical. We deform the tile T by moving to the right the left vertical 
boundary of T, until it hits the left vertical contour of a bubble, and vice 
versa for the right vertical boundary. We call T' and St' the new region 
and configuration obtained after the deformation. In passing from T, St to 
T', St' we increase the energy by at most 2|r|f'. Now we use the bound (2.37) 
which, as remarked after (2.37), is valid for all configurations consisting only 
of rectangular bubbles with the same orientation. The result is 


Et'^St') > edh*)\T'\ + CT^ + (2.42) 

n p n—1 

+^'^hi{edhi) -edh*)) + 7^^Wi{edwi) -edh*)), 

i=l ^ i=l 

where hi,... ,hn are the widths of the bubbles in St', and tci,..., Wn-i their 
separations. 

It T contains a hole, then either one of the h/s or tc/s is larger than £/5, 
or, the width of T' is smaller than 4i/5. In the first case, one of the terms 
ihi^edhi) — es{h*)) oi £wi{edwi)—es{h*)) is larger than (£^/5)|es(h*)|: to see 
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this recall that i > coh*, for a large enough constant cq, and use (2.30), which 


implies that es{h) is positive for h > i/5 in this case. In the second case, 
the difference between es{h*)\T'\ and es{h*)\T\ is larger than (£^/5)|es(h*)|. 
In both cases, we get a gain at least (f'^/5)|es(h*)|, which is larger than 
(const.To conclude the proof, note that under the stated 
assumptions on i (that is, coh* < i < (co|r|)“^ for a suitable constant cq), 
the error term (2 + C)|r|£ is smaller than c[nc{T) + , 

where c can be made as small as desired, by increasing cq (recall also that by 
definition of bad tile, either ndT) > 1/2, or XhoieiT) = 1). This concludes 
the proof. □ 


A The higher-dimensional case 

In this appendix we shall detail our main results in the case d > 3, and explain 
the main differences in their proof as compared to the two-dimensional case. 
The starting point is a representation of the energy in terms of droplets as 


in (2.3), whose boundaries, separating plus spins from minus spins, consist 


now of d — 1 dimensional plaquettes. Tiles are now d-dimensional cubes of 
side length i, and are used to divide space into good regions and bad cubes, 
with the good regions only containing “stripes” (i.e., quasi-one-dimensional 
regions of uniform spins, delimited by two flat parallel interfaces; they are 
slabs in d = 3), which can be oriented in d different directions. 


Our first claim concerns the fact that the localization bound (2.6) still 
holds, with the obvious notion of “corner”, namely d — 2 dimensional seg¬ 
ments where two plaquettes with different orientation meet. The proof of the 


analogue of (2.6) in higher dimensions is essentially the same as in d = 2, 
and relies on the analogue of (2.10), whose proof is in [201 App. D]. 


The key bound in Lemma for the good regions holds verbatim also for 
general d > 2. After the modihcation from G to G", each stripe has a dehnite 
width h but will not be a cuboid, in general; it is bounded by a union of d — 1 
dimensional cuboids Sk with width h and all other dimensions equal to i. As 
before, each Sk comes with two numbers, wi{sk) and W 2 {sk), measuring the 
distance to the next slice in the direction perpendicular to the stripes. The 


analogues of the slices gj introduced in Sect. 2.4.2 are cylinders with base 
area i'^~^ and various heights, which are obtained by adding up the various 
stripe widths hi and their separation Wi] they are oriented perpendicular to 
the stripes. With these modihed dehnitions, Eq. (2.17) still holds, with i 
replaced by i'^~^ in front of the hrst term on the right side, and f{wi, h,W 2 ) 
now denoting the interaction energy as depicted in Fig. with the upper 
strip of width h extended by i in the remaining d — 2 dimensions, while the 
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two lower ones are infinite in those directions. This function / satisfies the 
bound 


f{wi,h,W2) <i^ 


Co 


p—d—2 
111 - 

2 = 1,2 


(A.l) 


As already discussed in [201 App. A], the analogue of (|2.30) for general d is 




(A.2) 


from which it follows that es{h*) ~ |r|^^ and h* ~ |r| ^ for 

small r. Moreover, one easily deduces that (|A.l ) can be bounded by 


a 


wP 


< \t\w + (Cgirl ^)p-d-^w{es{w) - e^{h*)) , 


(A.3) 


which is the analogue of Eq. (2.34). The rest of the proof of Eq. (2.13) for 
general d remains unchanged. 

The analogue of Lemma for general d >2 takes the following form: 


Lemma 2’. For given d >2, there exist positive constants cq, C 2 and e such 
that, if —e < T < Q and c^h* <l< then the energy Et of any 

bad tile T eV satisfies 


Et{Bt) > i\{h*) + C2 [n,(T) + , (A.4) 


where Xhoie{T) is equal to 1 if T contains a hole, and 0 otherwise. 


Its proof is a rather straightforward adaptation of the one in d = 2, and 
we refer to [201 App. D], where the necessary changes were described in the 
case d = 3. 

Since every portion of the boundary of a good region Gi is adjacent to a 
bad tile, we have the bound ^ 2dC~^J\fB- In combination with 

the bounds above, this leads to the following generalization of Theorem 

Theorem 3’. For given d > 2, there exist positive constants Cq, Ci, e such 
that, if Jc — e < J < Jc and C^h* < i < {Cq{Jc — then for every 

sE {±1}^ and every finite set X C Z'^, 

Hxisxla*) > + + 

I J^G 

^ hj^h* i=l 

where Nc, and Gi are, respectively, the number of corners, the number 

of bad tiles containing a hole, and the good regions, associated with the infinite 
spin configuration a = {sxj^*x<=) coinciding with s^ on X and with a* on 
X^, defined via tiling with squares of side length i as described above. 


26 















Theorem implies the analogue of Theorem i.e., the fact that striped 
conhgurations with stripe width h* are inhnite volume ground states with 
trivial sectors for J close to Jc, and also the analogue of Theorem]^ stating 
that all infinite volume ground states that are invariant under translations 
by d — 1 lattice vectors are characterized by the existence of an interface 
separating the cubic lattice into two components, on each of which the 
conhguration is perfectly striped. 
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